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I.  STATEMENT  OF  THE  PROBLEM  STUDIED 


Since  the  Army  is  interested  in  techniques  of  producing  stronger  steels,  more  resistant  to 
fracture,  for  military  use  and  since  the  conventional  methods  of  testing  steels  for  strength  and 
fracture  toughness  are  destructive  in  nature,  it  was  proposed  that  this  project  be  a  fundamental 
study  to  understand  whether  mechanical  properties  can  be  correlated  with  magnetic  properties, 
and  if  so,  to  understand  how  and  why  this  happens  and  how  the  interrelationship  is  connected 
with  the  underlying  microstructure.  Because  harmonics  of  the  magnetic  induction  are  easily 
measured,  a  specific  interest  was  to  test  if  such  harmonics  could  be  used  to  nondestructively 
measure  tensile  strength  and  fracture  toughness  in  steel. 

Thus  the  project  was  (1)  to  investigate  the  interrelationship  between  mechanical  properties 
and  magnetic  properties  by  developing  an  understanding  of  how  microstructure  interconnects  the 
two  properties  and  (2)  to  test  if  nonlinear  harmonics  of  the  magnetic  induction  could  be  used  to 
nondestructively  determine  tensile  strength  and  fracture  toughness,  and  if  so,  how  and  why. 


II.  SUMMARY  OF  RESULTS 

A.  Effect  of  Tensile  Strength  on  Magnetic  Properties  and  on  Nonlinear 

Harmonic  (NLH)  Amplitudes 

Work  in  this  project  used  magnetic  hysteresis  modeling  to  explain  how  hysteretic 
magnetic  properties  are  correlatable  with  tensile  strength.  It  was  known  [1]  that  increasing  tensile 
strength  is  correlated  with  decreasing  grain  size  and  increasing  dislocation  density.  Certain 
parameters  in  magnetic  hysteresis  model  of  Sablik  and  Jiles  [2-4]  were  thus  modified  to  include 
the  contributions  of  grain  size  (average  grain  diameter)  and  dislocation  density.  Using  appropri¬ 
ate  functional  dependences  of  these  parameters  on  grain  size  and  dislocation  density,  a  computer 
study  was  then  done  to  investigate  how  these  attributes  affected  the  prediction  of  magnetic  prop¬ 
erties  like  coercivity,  remanence,  maximum  differential  penneability,  and  hysteresis  loss.  In  the 
range  of  grain  sizes  and  dislocation  densities  found  in  already  measured  specimens  [5],  the  com¬ 
puted  coercivity  was  found  to  show  a  linear  dependence  on  inverse  grain  size  and  a  propor¬ 
tionality  with  the  square  root  of  dislocation  density,  just  as  found  experimentally.  The  remanence 
and  maximum  differential  permeability,  on  the  other  hand,  decreased  with  respect  to  inverse 
grain  size  and  square  root  of  dislocation  density.  The  hysteresis  loss  showed  increase  with 
increasing  inverse  grain  size  and  square  root  of  dislocation  density,  in  agreement  with  experi¬ 
ment,  but  outside  the  range  tested  experimentally,  it  showed  opposite  tendencies  at  very  small 
grain  size  and  large  dislocation  density,  showing  the  effect  of  being  dominated  by  remanent 
behavior  instead  of  coercive  behavior.  These  results  were  reported  in  [5,6]. 

Modeling  previously  developed  [7]  for  computing  the  first  and  third  harmonic  ampli¬ 
tudes  was  modified  to  investigate  the  effect  of  grain  size  and  dislocation  density  on  the  first  and 
third  hannonic  amplitudes  of  the  magnetic  induction.  It  was  found  that  the  computed  harmonic 
amplitudes  decreased  with  decreasing  grain  size  and  increasing  dislocation  density.  Since 
decreasing  grain  size  and  increasing  dislocation  density  correlate  with  increasing  tensile  strength, 
in  effect  it  was  found  that  the  harmonic  amplitudes  decreased  with  increasing  tensile  strength,  as 
shown  experimentally  by  the  experimenters  at  IKPH  at  the  University  of  Hannover  in  Gennany. 
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Finally,  the  modeling  of  the  harmonic  amplitudes  was  extended  to  higher  frequen¬ 
cies,  taking  into  account  classical  eddy  current  power  loss  and  excess  power  loss  [8].  The  higher 
frequency  model  of  Jiles  [9]  was  modified  to  include  the  effects  of  grain  size  and  dislocation 
density.  A  numerical  solution  of  the  equations  had  to  be  employed.  The  results  [7]  agreed  with 
the  experimental  observation  that  at  a  fixed  frequency,  the  harmonic  amplitudes  decreased  with 
increasing  tensile  strength.  This  was  seen  for  specimens  with  coils  wrapped  around  the  speci¬ 
mens  themselves.  In  addition,  the  prediction  was  that  the  harmonic  amplitudes  were  smaller  at 
higher  frequencies  and  that  excess  loss  diminished  the  amplitudes  slightly  more.  The  decrease  of 
the  harmonic  amplitudes  with  increasing  frequency  was  seen  experimentally.  Experiments 
affirming  this  behavior  were  done  earlier  at  the  University  of  Hannover  [10]. 

In  this  project,  experiments  were  also  done  with  steel  sheets  and  air  coils  at  the  Uni¬ 
versity  of  Hannover  in  Gennany.  A  large  number  of  separate  specimen  sheets  were  tested  with 
respect  to  the  correlation  between  hannonic  amplitudes  and  tensile  strength.  The  tested  range  of 
tensile  strength  was  extended,  and  measurements  were  done  at  different  frequencies.  With  coils 
out  in  the  air  above  and  below  the  sheets,  an  unusual  feature  observed  was  that  the  third 
hannonic  amplitudes,  in  the  secondary  air  coils,  increased  instead  of  decreased  with  increasing 
tensile  strength.  This  differed  from  the  prediction  for  the  third  harmonic  amplitudes  found  inside 
a  specimen,  and  hence  this  effect  with  air  coil  configurations  had  to  be  further  investigated.  The 
effect  of  sheet  thickness  on  the  measurements  was  investigated.  It  was  found  that  amplitudes 
decreased  with  increasing  thickness  for  147  Hz,  and  for  higher  hannonics,  it  was  found  that  the 
amplitudes  decreased  more  sharply  and  then  leveled  off.  At  higher  frequencies  (>1400  Hz),  the 
hannonic  amplitudes  stayed  approximately  constant  with  increasing  thickness.  This  is  consistent 
with  the  expectation  that  at  higher  frequencies,  there  is  little  penetration  of  the  signal  into  the 
material,  so  the  amplitude  changes  little  as  thickness  is  further  increased.  At  lower  frequencies, 
the  signal  can  penetrate  the  sample  volume  much  more,  and  the  response  depends  on  the  volume 
percentage  of  sample  penetrated.  It  was  also  observed  that  tensile  strength  had  little  effect  on  the 
thickness  dependence. 

Meanwhile,  an  opportunity  developed  for  a  collaboration  with  Luc  Dupre  of  Ghent 
University  in  Ghent,  Belgium,  to  extend  our  predictions  of  microstructural  effects  on  magnetiza¬ 
tion  to  the  predictions  of  the  microstructural  dependence  of  Preisach  parameters  used  to  fit 
hysteresis  loops  by  what  is  known  as  Preisach  hysteresis  modeling.  Usually,  the  Preisach  model 
is  used  to  predict  hysteresis  loops  based  on  fitting  sequences  of  experimental  inner  hysteresis 
loops  until  the  final  experimental  saturation  hysteresis  loop  is  reached.  In  steels,  as  part  of  the 
fitting  process,  a  Lorentzian  Preisach  distribution  function  is  used  to  do  this.  Dupre  had  studied 
the  Lorentzian  Preisach  distribution  functions  needed  to  fit  the  experimental  hysteresis  curves  of 
a  family  of  steel  specimens  of  different  grain  sizes  and  had  determined  how  the  various  Lorent¬ 
zian  Preisach  distribution  parameters  varied  with  grain  size  [11].  In  a  new  approach  to  the 
modeling,  hysteresis  curves  of  inner  and  saturation  loops  were  obtained  using  the  Jiles- Atherton- 
like  [12]  approach  that  had  been  used  in  this  project.  Hysteresis  loops  were  computed  first  for  a 
set  of  different  grain  sizes  and  constant  dislocation  density,  and  then  for  a  set  of  different  dislo¬ 
cation  densities  and  constant  grain  size.  The  first  set  of  loops,  involving  different  grain  sizes,  was 
fitted  using  a  Lorentzian  Preisach  distribution  function.  Results  exhibited  good  agreement  with 
the  earlier  fit  by  Dupre  [11]  to  experimentally  obtained  loops  for  different  grain  sizes.  The  other 
set  of  computed  loops,  for  varying  dislocation  densities,  was  used  to  predict  how  the  Lorentzian 
Preisach  distribution  parameters  would  behave  for  a  new  set  of  experimental  specimens  with 
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varying  dislocation  density  and  constant  grain  size.  These  calculations  are  now  available  for 
future  experimental  tests.  Physical  arguments  given  in  our  published  paper  [13]  offer  explana¬ 
tions  as  to  why  the  Preisach  parameters  behave  the  way  they  do. 

Through  conference  discussion,  it  was  found  that  there  existed  unpublished  data  for 
hysteresis  loss  conflicting  with  predictions  of  our  model  in  the  case  of  very  small  grain  size  and 
large  dislocation  density.  F.  J.  G.  Landgraf,  in  the  same  conference  discussion,  suggested  what 
might  be  responsible  for  the  apparent  conflict.  He  noted  that  hysteresis  loss  experimenters  com¬ 
pare  hysteresis  loops  of  the  same  maximum  flux  density  Bmax  whereas  in  our  modeling, 
hysteresis  loops  of  the  same  maximum  magnetic  field  Hmax  were  being  compared.  The  computer 
program  was  thus  restructured  to  output  hysteresis  loops  all  with  the  same  maximum  flux 
density.  This  time,  predicted  hysteresis  loss  revealed  a  linear  increase  with  increasing  inverse 
grain  size  and  increasing  square  root  of  dislocation  density,  just  as  found  experimentally  for  very 
small  grain  and  large  dislocation  density.  So,  the  model  does  indeed  reproduce  the  experimental 
dependencies.  Remanent  field  behavior  in  this  case  does  not  influence  hysteresis  loss  away  from 
the  pattern  followed  by  coercive  field  behavior  because  the  remanent  field  remains  essentially 
constant  for  all  loops  of  the  same  maximum  flux  density.  Two  papers  were  prepared  discussing 
these  effects  [14,15].  Also  discussed  in  these  papers  is  the  effect  of  uniaxial  microcrystalline 
anisotropy  on  the  microstructural  dependence  of  the  magnetic  properties.  One  byproduct  of  these 
new  papers  is  a  mathematical  argument  justifying  microstructural  interrelationships  among  the 
Jiles- Atherton  hysteresis  parameters  found  earlier. 

Earlier  we  had  discussed  results  indicating  that  with  primary  and  secondary  coils 
wrapped  around  the  specimen,  the  nonlinear  hannonic  (NLH)  amplitudes  of  the  magnetic  induc¬ 
tion  decreased  with  increasing  tensile  strength.  We  have  now  determined  by  finite  element 
modeling  (See  Section  C)  that  the  configuration  to  be  used  in  a  steel  plant  shows  somewhat 
different  behavior.  For  the  steel  plant,  air  coils  (primary  and  secondary)  are  to  be  used  above  and 
below  a  moving  steel  sheet.  The  secondary  coil  voltage,  however,  is  a  measure  of  the  flux  den¬ 
sity  in  the  air,  and  not  in  the  specimen  itself.  The  flux  density  in  the  air  is  modified  by  the 
presence  of  the  steel  sheet,  but  in  the  secondary  air  coils,  hannonics  change  somewhat  differently 
from  the  flux  density  harmonics  that  would  be  measured  in  the  specimen  itself.  For  the  air  coils, 
the  third  harmonic  amplitude  monotonically  increases  with  increasing  tensile  strength,  and  other 
hannonics  show  mixed  behavior,  with  initial  decrease  and  then  increase.  A  way  to  look  at  this  is 
to  realize  that  with  flux  concentrated  in  the  specimen,  it  gets  more  difficult  to  change  the  flux  in 
the  specimen  with  increasing  tensile  strength  and  increasing  coercivity,  and  thus  the  air  coils 
tend  to  drive  mostly  flux  changes  in  the  air,  which  are  more  easily  driven,  and  hence  the  higher 
harmonics  increase  at  the  higher  tensile  strengths. 

B.  Effect  of  Fracture  Toughness  on  Nonlinear  Harmonic  Amplitudes 

Because  of  the  retirement  of  co-P.I.  Dieter  Stegemann,  work  at  Hannover  was  con¬ 
cluded  in  June  2001,  and  it  fell  to  new  co-P.I.  Gary  Burkhardt  at  SwRI,  who  is  experienced  at 
nonlinear  harmonic  (NLH)  nondestructive  evaluation  (NDE)  [16,17],  to  continue  experimental 
work,  this  time  investigating  the  effect  of  changing  fracture  toughness  on  the  harmonic 
amplitudes. 


Experimentally,  it  was  decided  that  the  NLH  measurements  would  be  conducted  on 
Charpy  specimens  since  the  Charpy  test  is  one  way  of  measuring  fracture  toughness.  Fracture 
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toughness  was  varied  by  heat-treating  the  specimens  prior  to  testing,  according  to  standard 
procedures.  In  setting  up  such  procedures,  we  have  been  guided  by  one  of  the  metallurgists  at 
SwRI,  N.  S.  Cheruvu.  Experiments  were  carried  out  on  two  different  types  of  steel,  namely  4340 
and  9-4-30  steel.  Once  heat  treatments  were  carried  out,  the  specimens  were  machined  with  V- 
notches,  according  to  the  standard  Charpy  dimensional  specifications  [18],  and  then  NLH 
measurements  were  performed  at  two  separate  frequencies  for  each  specimen.  Finally,  the 
destructive  Charpy  tests  were  carried  out  and  yield  strength  was  measured.  In  doing  the  standard 
Charpy  measurements,  a  knife  edge  pendulum  was  dropped  from  a  height  h.  The  minimum 
height  from  which  the  specimen  fractures  gives  a  measure  of  the  fracture  energy,  nigh.  Speci¬ 
mens  fractured  at  different  temperatures  fracture  with  different  fracture  energies.  The  fracture 
energy  vs.  temperature  curve  has  a  step-like  shape,  going  from  a  lower  brittle  fracture  energy  to  a 
higher  ductile  fracture  energy,  known  as  the  upper  shelf  energy  CVN.  Different  heat  treatments 
shift  the  mean  temperature  at  which  the  step  occurs,  and  also  different  CVN  obtain  for  different 
specimens.  A  relationship  exists  between  CVN,  the  yield  strength  oy,  and  stress  intensity  factor 
Kjc,  which  is  often  taken  as  a  measure  of  fracture  toughness.  The  relationship  is  the  following 
[19]: 

CA/Kf  =(5A){cw-K/20)}  (d 

In  this  relationship,  ay  is  in  ksi,  CVN  is  in  ft-lb,  and  K\c  is  in  ksi  V in.  We  have  used  this 
relationship  to  compute  K\c.  However  the  Charpy  curves  for  fracture  energy  vs.  temperature  for 
each  specimen  often  had  different  temperatures  for  the  beginning  of  the  upper  shelf  compared  to 
the  upper  shelf  of  the  curve  for  %  ductile  fracture  vs.  temperature.  Thus  it  was  not  clear  as  to 
where  to  position  the  upper  shelf  and  how  to  compute  the  CVN.  We  took  a  simple  estimate  and 
designated  the  fracture  energy  at  room  temperature  to  be  the  CVN.  Clearly,  there  are  errors  in 
such  an  estimate.  We  present  two  tables.  In  Table  1,  we  show  the  tempering  temperature,  target 
hardness,  measured  hardness,  and  range  of  K\c  expected  for  4340  steel  specimens  at  the  indicated 
tempering  temperatures.  In  Table  2,  we  show  for  each  4340  and  9-4-30  specimen  the  estimated 
CVN,  the  measured  oy,  and  the  computed  Klc  (averaged  over  two  specimens  for  each  specimen 
type).  Figures  1  and  2  display  the  dependence  of  the  third  and  fifth  hannonics  for  each  K\c  for 
the  two  different  specimen  types,  with  150-Hz  signal  applying  to  Figure  1  and  1500  Hz  applying 
to  Figure  2.  Clearly,  our  results  do  not  indicate  a  correlative  relationship  between  NLH  ampli¬ 
tude  and  Kic  for  150  Hz,  and  indicate  only  at  most  a  weak  correlation  between  NLH  amplitude 
and  Kic  for  1500  Hz.  Since  the  K\c  values  that  we  have  used  are  probably  in  error,  it  is  clear  that 
a  more  acceptable  approach  would  be  to  measure  K]c  directly  by  a  more  standard  measuring 
technique  [20], 

In  examining  the  microstructural  features  that  contribute  to  fracture  toughness,  a 
good  correlation  between  fracture  toughness  and  microstructural  features  was  not  obtainable  also 
with  the  modeling,  at  least  as  yet.  It  so  far  appears  that  different  features  contribute  conflicting 
effects.  So  the  inconclusive  correlation  obtained  experimentally  between  fracture  toughness  and 
NLH  amplitude  may  exist  even  with  better  measurement  of  K\c.  Our  retired  co-PI  has  recently 
also  done  some  measurements  to  explore  potential  correlation  between  NLH  and  fracture  tough¬ 
ness  properties,  and  he  too  has  found  inconclusive  correlation  [21].  Thus,  it  may  well  be  that, 
although  NLH  gives  good  correlation  between  tensile  strength  and  NLH  amplitude  (even  in  the 
case  of  air  coils,  where  trends  differ  from  that  obtained  with  wrap-around  coils),  NLH  may  not 
be  a  good  tool  for  NDE  evaluation  of  fracture  toughness.  More  work  is  needed  to  establish  this. 
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Table  1.  Information  about  Sample  Preparation  for  4340  Steel 

4340  STEEL 


Table  2.  Computation  of  K\c  for  Both  Steels 
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Figure  1.  Third  and  fifth  harmonic  amplitudes  vs.  Kic  for  4340  steel 
and  9-4-30  steel  for  150-Hz  signal 
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Figure  2.  Third  and  fifth  harmonic  amplitudes  vs.  K\c  for  4340  steel 
and  9-4-30  steel  for  1500-Hz  signal 
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C.  Technology  Transfer 

As  mentioned  in  our  interim  reports,  the  technical  group  at  the  University  of  Hannover  has 
had  additional  support  from  the  German  steel  industry  (namely,  Salzgitter  AG)  to  develop  pilot 
instrumentation  to  be  used  on-line  in  the  factory  to  test  sheet  steel  for  its  mechanical  properties — 
in  particular,  tensile  strength.  An  air  coil  configuration  has  been  built  for  this  purpose  and  is 
being  tested  at  the  steel  plant. 

A  separate  European  Union  subcontract  was  issued  to  SwRI  to  use  finite  element  modeling 
and  the  magnetic  model  developed  for  tensile  strength  effects  to  unravel  the  physics  of  what 
happens  with  a  complicated  coil  configuration  and  to  help  optimize  the  coil  design.  Thus  the 
present  ARO  project  has  resulted  in  additional  work,  which  should  bring  about  better  technology 
transfer. 

Another  application  has  been  to  develop  a  magnetic  test  for  determining  whether  a  weld  is 
annealed  [22].  The  modeling  on  grain  size  effects  and  dislocation  density  effects  on  magnet¬ 
ization  was  used  as  input  into  a  finite  element  calculation  predicting  magnetic  signal  changes  as 
a  probe  is  moved  across  the  various  regions  of  a  weld — base  metal  (away  from  the  weld),  heat- 
affected  zone  (which  exhibits  increasing  grain  size  and  increasing  dislocation  density  as  one 
moves  toward  the  weld  center),  and  finally  the  fusion  zone  in  the  central  part  of  the  weld  with 
large  grain  sizes,  but  also  large  dislocation  densities  if  the  weld  is  unannealed.  It  turns  out  that 
the  dislocation  density  effect  dominates  the  grain  size  effect  when  the  weld  is  unannealed,  and 
the  permeability  (determined  from  the  magnetic  signal)  decreases  sizably  as  one  moves  from  the 
base  metal  to  the  central  part  of  the  weld.  On  the  other  hand,  when  the  weld  is  annealed,  disloca¬ 
tion  densities  greatly  decrease,  and  thus  the  grain  size  effect  dominates  the  dislocation  density 
effect,  and  the  permeability  actually  increases  slightly  or  moderately  as  one  moves  from  the  base 
metal  to  the  central  part  of  the  weld.  The  two  different  behaviors  tell  you  immediately  whether 
the  weld  is  annealed  or  not.  This  then  is  an  example  of  how  the  basic  physics  exposed  by  this 
project  can  be  utilized  to  design  an  eminently  usable  NDE  test  involving  weld  characteristics  in 
steels,  and  thus  points  to  yet  another  technology  transfer  application. 
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Abstract—  Microstructural  attributes  such  as  grain  size  d  and 
dislocation  density  Q  affect  the  hysteretic  magnetic  properties  of 
steels  because  they  affect  domain  wall  movement  and  pinning.  In 
an  earlier  paper,  a  model  was  proposed  for  computing  hysteresis 
loops  based  on  the  effect  of  grain  size  and  dislocation  density.  In 
that  paper,  hysteresis  loops  were  compared  that  all  had  the  same 
maximum  field  Hmax.  The  result  was  that  coercivity  departed 
from  linear  relationships  with  inverse  grain  size  (viz.  lid)  and 
Cd'1  for  large  values  of  lid  and  £/1/2-  The  same  was  true  of 
hysteresis  loss  WH,  except  that  hysteresis  loss  even  showed  a  peak, 
first  increasing  and  then  decreasing  with  increasing  ltd  and  £//2. 
This  kind  of  behavior  had  not  been  seen  by  experimenters, 
particularly  core  loss  people.  It  was  learned  that  the  core  loss 
experimenters  compared  hysteresis  loops  of  the  same  maximum 
flux  density  Bmax  instead  of  the  same  Hmax.  In  this  paper,  we  use 
the  model  previously  formulated  to  produce  hysteresis  loops  with 
the  same  Bmax.  Indeed,  the  appropriate  linear  relationships  are 
found.  The  paper  also  addresses  effects  of  uniaxial  anisotropy  on 
these  microstructural  magnetic  effects  and  why  two  hysteresis 
parameters  are  affected  by  microstructural  variation. 

Index  Terms — Hysteresis  modeling,  magnetic  materials,  micro- 
structural  effects,  uniaxial  anisotropy  effects. 

Introduction 

IN  steel,  two  important  microstructural  features  that  affect 
magnetic  hysteresis  are  (1)  grain  size  (i.e.,  average  grain 
diameter  d)  and  (2)  dislocation  density  Q.  Since  domain  walls 
tend  to  pin  at  grain  boundaries,  the  pinning  of  domain  wall 
motion  increases  with  increasing  total  grain  boundary  length 
as  grain  size  d  decreases.  Since  coercivity  Hc  reflects  amount 
and  strength  of  pinning,  we  expect  Hc  to  increase  as  d 
decreases.  Similarly,  as  dislocation  density  Q  increases, 
dislocations  begin  to  get  entangled,  forming  strong  pinning 
centers  for  domain  walls,  so  impeding  domain  wall  motion. 
Thus,  as  Q  increases,  so  also  does  Hc.  It  has  been  generally 
established  experimentally  [1,2]  that  Hc  has  a  linear 
relationship  with  respect  to  A  +  Bid  and  <£//2,  where  A  and  B 
are  constants.  A  similar  relationship  was  also  generally  found 
for  hysteresis  loss  WH  [3,4]. 

Despite  all  this  experimental  work,  no  one  had  actively 
modeled  the  effect  of  grain  size  and  dislocation  density  on 
magnetic  hysteresis.  In  two  recent  papers  [5,6],  this  changed 
when  one  of  the  present  authors  (M.S.)  presented  a  hysteresis 
model  which  was  a  modification  of  an  earlier  hysteresis  model 
due  to  Jiles  and  Atherton  [7]  The  model  successfully  exhibited 
linear  behavior  with  A  +  B/d  and  Q1/2  for  grain  sizes  of  order 
larger  than  15/v  (  where  //  represents  1  micron  =  10"6  m  )  and 
for  smaller  dislocation  densities.  In  comparing  hysteresis 


loops,  the  loops  were  all  taken  to  the  same  maximum  magnetic 
field  Hmax.  The  result  appeared  to  agree  with  experimental 
results,  since  all  the  published  experiments  were  restricted  to 
grain  sizes  larger  than  15//.  Core  loss  experimenters  have  now 
privately  indicated  that  they  can  find  linear  results  for  smaller 
grain  sizes  than  15//.  However,  they  compare  hysteresis  loops 
taken  to  a  constant  Bmax  instead  of  a  constant  Hmax. 

It  was  decided  therefore,  in  this  paper,  to  compare  modeled 
hysteresis  loops,  all  taken  to  the  same  Bmax.  It  was  reasoned 
that  for  constant  Bmax,  the  remanence  B,  would  remain 
approximately  constant  regardless  of  grain  size  and  dislocation 
density.  Thus,  hysteresis  loss  WH,  which  is  approximately 
B,HC,  should  show  the  same  pattern  of  behavior  as  Hc  On  the 
other  hand,  in  the  case  where  all  loops  are  taken  to  a  constant 
Hmax,  the  remanence  decreases  with  increasing  inverse  grain 
size  and  dislocation  density,  whereas  the  coercivity  increases, 
causing  a  competition  between  B,  and  Hc  in  contribution  to 
Wh,  with  a  maximum  resulting  in  Wu  due  to  the  competition. 
We  thus  expect  much  more  regular  behavior  for  loops  of 
constant  Bmax. 

The  Model 

We  start  by  referring  to  the  basic  hysteresis  model  of  Jiles 
and  Atherton  [7],  which  has  been  modified  to  include  the 
effect  of  stress  by  Sablik  and  Jiles  [8]. 

In  the  Jiles-Atherton  model,  the  total  magnetization  M  is  the 
sum  of  a  reversible  (Mrev)  and  an  irreversible  (Mirr)  compo¬ 
nent.  These  components  are  given  by 

Mrev=c  (Ma  -  Mjrr )  (1) 


Mo  dHe 

Here,  Ma  is  the  anhysteretic  magnetization,  given  as 

Ma(He)  =  MsL(He/a),  (3) 

where  L(x )  =  coth  x  -  1/x  is  the  so-called  Langevin  function, 
and  where  He  is  the  effective  magnetic  field  in  the  material, 
i.e. 

He=  H  +  aMa  .  (4) 

The  five  parameters  Ms,  c,  a,  k  and  a  are  all  parameters  of  the 
material.  The  parameter  5  is  +1  or  -1,  depending  on  whether 
H  is  increasing  or  decreasing.  Equation  (2)  can  be  reexpressed 
as  a  differential  equation  for  dMjrrldH  [7,8]. 
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Since  k  mathematically  controls  the  amount  of  hysteresis 
that  is  present,  it  is  proportional  to  the  coercivity  and  hence 
has  the  same  dependences  as  the  coercivity.  Thus,  we  write 

k  =  \Gx  +  G2/d\  i/f  K-  (5) 

For  d  =  d*  =  20//  and  Q  £,*  —  1  x  1012  /m2,  we  choose  G\ 
and  Gi  so  that 

(< Gl+Gjd' )  {Cf  =  1.  (6) 

A  choice  that  satisfies  this  is  G2  =  10  x  10"12  m2  and  Gt  =  0.5  x 
10"6  m.  Other  choices  for  G2  and  G\  also  satisfy  (6).  We  find 
that  various  choices  for  G2  and  G\  correlate  with  the  amount 
of  spread  in  the  hysteresis  behavior  due  to  variation  of  grain 
size  and  dislocation  density.  Note  that  when  d  and  Q  satisfy 
(6)  (i.e.,  when  d  =  d*  and  Q  =  Q*),  then  k  =  k0.  The  20 // 
choice  for  grain  size  d*  represents  a  typical  grain  size  that  has 
been  investigated.  The  choice  of  Q*  =  1  x  1012/m2  represents  a 
dislocation  density  quoted  by  papers  on  plastic  deformation  as 
a  typical  value  for  dislocation  density  in  undeformed  steel 
[4,9].  Other  papers  [1,10]  have  quoted  values  of  the  order 
10’%;2,  and  such  values  were  used  in  [3]  and  [4].  Because  the 
range  of  Q  used  here  is  now  of  order  10I2/to2,  the  values  of  G\ 
and  G2  are  altered  to  accommodate  this  range. 

Scaling  constant  a  is  proportional  to  domain  density  in  the 
demagnetized  state  [8],  which  is  determined  by  pinning  site 
density,  in  turn  proportional  [7]  to  pinning  constant  k.  Thus,  a 
has  the  same  dependence  on  d  and  Q  that  k  has,  and  hence 

a  =  {[G3+Gjd]cf}  ao.  (7) 

We  shall  show  mathematically  that  a  is  directly  related  to  Hc. 
If  we  define  G3  and  G4  in  the  same  way  as  Gy  and  G2,  using 
(6)  for  d  =  d  and  Q  =  Q  ,  then  we  can  set  G3  =  Gj  and  G4  = 
G2. 

Results 

In  the  following  analysis,  we  have  restricted  the  grain  size  d 
to  10//,  15 //,  20//,  25//,  and  30//.  Also,  we  have  restricted  the 
dislocation  density  Q  to  0.25,  0.49,  1,  2.25  and  4  xl012/m2. 
The  two  most  extreme  conditions  occur  for  microstructural 
value  sets  of  (4  x  1012/to2,  10// )  and  (0.25  x  1012  /m2,  30// ). 
For  hysteresis  loops  all  of  the  same  Bmax,  the  first  value  set 
corresponds  to  the  loop  with  largest  Hmax,  and  the  second 
value  set  corresponds  to  the  loop  with  smallest  Hmax.  By 
keeping  Bmax  constant  for  all  the  hysteresis  loops,  we  end  up 
with  loops  of  widely  varying  shape  and  breadth,  as  we  vary  d 
and  Q.  To  arrive  at  Bmax,  one  increments  FI  in  smaller  and 
smaller  increments  as  Bmax  is  neared.  Flere,  we  have  Bmax=  1.03 
T.,  c  =  0.25,  ko//Jo=  1200  A/m  (where  /u0  is  the  permeability  of 
free  space),  a0=  1100  A/m,  a=  8.44  x  10~6  (where  a  is  related 
to  As  via  (30)  in  [8]),  and  Ms  =  1.585  x  106  A/m.  Also,  we  have 
used  the  values  for  G2  and  Gj  given  just  below  (6).  The 
parameters  are  thus  chosen  for  hard  magnetic  material.  From 
the  various  hysteresis  loops,  we  obtain  the  magnetic  properties 
of  coercive  field  //,,  relative  permeability  lu'li,,  at  field  H  =  Hc, 


remanent  flux  density  Br  and  hysteresis  loss  Wrl  for  each  loop 
in  the  set  of  25  loops  corresponding  to  the  different  (Q,  d). 

Fig.  1  is  a  set  of  plots  of  the  different  magnetic  properties 
against  square  root  of  the  dislocation  density.  Two  things  are 
striking.  The  coercive  field  Hc  is  directly  proportional  to  the 
square  root  of  the  dislocation  density  for  all  values  of  Qm, 
exhibiting  different  slopes  for  different  grain  sizes.  There  is  no 
deviation  from  linear  proportionality,  as  is  seen  for  the  curves 
when  loops  of  the  same  Hmax  [5]  are  compared.  Thus,  the 
better  way  to  compare  hysteresis  loops  for  different  grain  sizes 
and  dislocation  density  is  to  compare  loops  that  are  all  taken  to 
the  same  Bmax.  The  same  is  true  for  hysteresis  loss  Wh .  The 
hysteresis  loss  is  also  proportional  to  the  square  root  of  the 
dislocation  density  in  this  model  for  all  values  of  i/d12,  with 
different  slopes  for  different  grain  sizes.  No  maximum  in  Wh, 
appears,  as  was  found  when  loops  of  the  same  Hmax  were  com¬ 
pared  [5],  The  relative  permeability  at  Hc  does  not  exhibit 
linear  behavior  and  decreases  nonlinearly  with  increasing  i/d12. 
This  is  similar  to  what  was  observed  [5]  for  loops  all  with  the 
same  Hmax.  The  remanent  Br  is  approximately  the  same  for  all 
values  of  i/d  and  d.  as  was  anticipated.  (See  the  Introduction.) 

We  also  look  at  plots  of  the  magnetic  properties  against 
inverse  grain  size  1/d,  as  shown  in  Fig.  2.  Flere,  Hc  and  WH 
vary  essentially  linearly  with  inverse  grain  sizes,  but  with 
plots  for  different  Q  having  different  intercepts.  The  relative 
permeability  at  Hc  varies  nonlinearly,  decreasing  with 
increasing  \/d.  The  remanent  flux  density  is  approximately 
constant. 

Inclusion  of  Uniaxial  Anisotropy 

One  can  modify  the  basic  model,  as  outlined  in  ( 1)— (4),  to 
include  uniaxial  anisotropy.  The  key  modification  is  that  [11] 

He=H  +  (a-Ku)Ma,  (8) 

where 

yu=2^/(//0M2),  (9) 

where  Ku  is  a  measure  of  the  uniaxial  anisotropy  energy  Eu  per 
unit  volume.  Equations  (l)-(3)  still  apply. 

In  this  section,  we  study  effects  of  uniaxial  anisotropy,  due 
perhaps  to  grain  orientation,  on  magnetic  hysteresis  properties. 
For  this  case,  Ms  =  1.585  x  106  A/m,  ko/u0  =  500  A/m,  a0=  600 
A/m,  c  =  0.25  and  Gj  =  0.81  x  10~6  and  G2  =  3.8  x  10~12  to2. 
We  treat  the  cases  of  Ku  =  -1200  //to3,  0,  and  +1200  //to3, 
with  negative  Ku  corresponding  to  uniaxial  anisotropy  and 
positive  Ku  to  perpendicular  planar  anisotropy.  The  interesting 
plots  are  those  of  magnetic  properties  against  i/d'2 ■  Fig-  3 
shows  these  plots  for  Hc  and  Wh .  In  particular,  while  the  plots 
of  Hc  vs.  i/d'1  and  Wh  vs.  t/d /2  still  exhibit  linear  behavior,  the 
plots  for  the  different  grain  sizes  extrapolate  to  a  value  along 
the  ordinate  axis  that  is  nonzero  when  Ku  is  nonzero.  In  par¬ 
ticular,  when  Ku  =  +1200  //to3  (and  with  Bma=  1.62  T  for  the 
loops),  the  value  to  which  the  Hc  and  WH  plots  extrapolate  is 
positive  and  nonzero,  and  so  the  Hc  and  Wh  cannot  ever  be 
zero  in  this  anisotropic  case.  K„=0  lines  all  have  intercept 
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Hc=0  and  Wh=0.  On  the  other  hand,  when  Ku  =  -1200  J/m3 
(and  with  Bmax=lJ0  T  for  the  loops),  the  intercept  to  which  Hc 
and  WH  extrapolate  is  negative,  which  actually  cannot  be  a 
physical  value  because  Hc  and  WH  can  only  be  positive.  The 
extrapolated  lines  also  intercept  the  abscissa  axis.  All  values 
of  the  abscissa  £//2  less  than  the  extrapolated  value  at  the 
abscissa  intercept  are  values  for  which  the  hysteresis  model  is 
unstable  and  for  which  it  cannot  be  used.  Since  uniaxial 
anisotropy  can  lead  to  discontinuous  transitions,  it  is  not 
surprising  that  continuous  hysteresis  will  have  some 
unrealizable  ranges. 

To  see  why  the  ordinate  intercept  is  nonzero  when  Ku  is 
nonzero,  remember  that  the  Langevin  function  L(x)  =  coth(x) 
-  1/x,  and  recall  that  the  series  for  coth(x)  is  such  that 

L(x)«x/ 3 -x3/45  + 2*7945  «x/3,  (10) 

and  since  x  =  H/a,  we  have  that 

Ma/Ms«[H  +  (a-Kll)Ma]/3a .  (11) 

If  H=HC,  then  from  (11),  it  follows  that 

Hc  *  3  aMa  (7/  )/Ms -(a-  k„  )Ma  (Hc )  (12) 

If  Ku  is  negative,  the  anisotropy  term  ku  Ma(Hc)  subtracts 
from  the  Ku=  0  value  of//,.  lfKu  >  0,  it  adds.  Also  note  that  if 
pinning  constant  ko  tends  to  zero,  then  M  tends  to  Ma  and  since 
M  =  0  at  H  =  Hc,  it  follows  that  with  Ku  =  0,  then 

Hc  *  ((3a/ AT  )-a)Ma{Hc)»  ((3  a/Ms  )-a)M(Hc)  =  0.(13) 

Thus,  the  Hc  intercept  is  zero  for  K„  =  0;  nonzero  otherwise. 

If  one  uses  (7)  for  a ,  then 

V  *  {(3a„ /M,  )[G1+G2 /d]  CJ2  ~(a~  k„  )} Ma  (Hc ) .  ( 1 4) 

Thus,  consistent  with  model  and  experiment,  it  clearly  is 
seen  that  Hc  is  linearly  dependent  on  lJ  1  and  on  [Gi  +  G2/d  ]. 
This  behavior  derives  from  a,  and  hence  a  must  behave  like  k 
does,  as  asserted  earlier,  using  a  different  argument. 
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Fig.  1.  Magnetic  properties  against  square  root  of  dislocation  density  Q,  with 
material  parameters  set  as  in  the  text.  Shown  are  coercive  field  Hc,  relative 
permeability  ju/jUq  at  Hc,  remanent  flux  density  Br  and  hysteresis  loss,  WH. 
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Fig.  2.  Plots  of  magnetic  properties  against  inverse  grain  size.  The  plots  for 
each  value  of  dislocation  density  correspond  to  straight  lines  with  different 
intercepts,  in  the  case  of  Hc  and  Wh .  Here,  parameter  values  are  as  in  Fig.  1 . 
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Fig.  3.  Plots  of  (a)  coercive  field  and  (b)  hysteresis  loss  against  C,^12  for 
different  grain  sizes  and  different  anisotropies  Ku.  The  text  shows  why  straight 
lines  for  the  different  anisotropies  converge  to  different  intercepts. 


